Tunnelling matrix elements with antiferromagnetic Gutzwiller wave functions 
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We use a generalized Gutzwiller Approximation (GA) elaborated to evaluate matrix elements 
with partially projected wave functions and formerly applied to homogeneous systems.— 1 ^ In the 
present paper we consider projected single-particle (hole) excitations for electronic systems with 
antiferromagnetic (AFM) order and obtain the corresponding tunnelling probabilities. The accuracy 
and the reliability of our analytical approximation is tested using the Variational Monte Carlo 
(VMC). Possible comparisons with experimental results are also discussed. 

PACS numbers: 



I. INTRODUCTION 



Strongly correlated electronic systems can be stud- 
ied by means of several methods like quantum Monte 
Carlo^— , dynamical mean-field theory^, numerical and 
functional renormalization group£~— , density functional 
theorj*i2r— , Gutzwiller method 13 and slave bosonsii 4 - 
Most approaches for the strong coupling regime are nu- 
merical: they can explore cases difficult to be treated 
analytically but they are typically computer-time con- 
suming and not always allow a full physical insight into 
the investigated quantities. On the other side, analyti- 
cal computations lead to results in form of explicit re- 
lations but they often require some approximations or 
ansatz not totally under control in their accuracy. This 
is also the case of the variational studies which use the 
Guztwiller wave functions and evaluate the observables 
via the Gutzwiller approximation. We focus our paper 
on the evaluations of matrix elements with Gutzwiller 
functions, motivated by the past success of studies with 
these projected wave functions. In fact, they provided 
significant progresses in the characterization of strongly 
correlated systems^ 5 ., also to study the effective models 
proposed to explain the properties of high-temperature 
superconductors We also intend to establish more 
precisely the reliability and the drawbacks of our approxi- 
mated analytical computations: to this aim, we calculate 
some test quantities also with the VMC, which provides 
exact results within the statistical uncertainties. 

Gutzwiller proposed to study electronic systems with 
repulsive onsite interactions using projected wave func- 
tions of the form \^f) = P\ty ), where P is a projector 
which suppresses the double occupancies of the uncor- 
rected trial state |\J/ )iiii If the latter is the homoge- 
neous state, the projector is P = Yl^l ~ n it n il) where 
rii a = c icr Ci a , with c\ a (ci a ) the creation (annihilation) 
operator for an electron of spin a on a site i. 

In the GA, one evaluates the expectation values of op- 
erators O using the projected wave functions ^. The 
effects of the projection on the state \^o) are approx- 
imated by a classical statistical weight factor g a which 



multiplies the unprojected quantum result; 13 ' 15 
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Basically, the projection operator P reduces the num- 
ber of the allowed states in the Hilbert space, and in the 
simplest approximation this renormalization is taken into 
account through combinatorial factors. In order to ob- 
tain the renormalization factors, one counts the number 
of states that contribute to (*|6|*) and to (^oM^o), 
and their ratio provides the factor. If for example O is 
the kinetic operator TL^ = cj c. 
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sponding factor g tcr expresses the reduction of the kinetic 
energy for correlated electrons. Such renormalization fac- 
tors are obtained for example in Ref . Il3lll5lll9l - l23l and also 
in equivalent slave-boson approaches. 14,24 In our present 
work, we use Gutzwiller factors obtained neglecting any 
noncombinatorial configuration dependences and inter- 
site correlations. For further discussions, the interested 
reader can consider Refs. [l9ll20l for example. As we will 
explain better in Sect. II, we also introduce an auxiliary 
Gutzwiller projector P' to handle computations with par- 
tially projected wave functions: P' suppresses the dou- 
ble occupancies for all the lattice sites but for a reservoir 
sitej 1 ^ 

In this paper we consider projected single-particle 
(hole) excitations and evaluate the corresponding tun- 
nelling probabilities for an AFM, i.e. a system on a 
bipartite lattice with collinear spin order and staggered 
magnetization m. We consider the AFM case, also be- 
cause this is a recurring wave function for the ground 
state of (lowly-doped) strongly correlated systems. 

We emphasize that in our approach the results are due 
only to the physical properties of the trial state and not 
to the choice of a specifical Hamiltonian: in this sense, 
they are model-independent but they are not universal, 
because they rely on the features of the chosen Gutzwiller 
wave function (Fermi sea, BCS superconductor, AFM...) 

Our results can be helpful to understand the spectra 
probed in one-particle spectroscopy experiments (tun- 
nelling, photoemission) for materials displaying a long- 
range collinear AFM order. Actually, to date there are 
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not many tunnelling experiments dedicated for strongly 
correlated materials in their long-range AFM region. 
In fact, the standard Scansion Tunnelling Microscopy 
(STM) method can not detect spin structures 2 ^ and the 
improved Spin-Polarized STM can not be applied for in- 
sulating materials i 26 ' 27 In order to study the AFM or- 
der in correlated insulators, a new spin-probe method 
has been developed: the Magnetic Exchange Force 
Microscop y 26 ' 28 , which does not involve one-particle pro- 
cesses. We will discuss possible future comparisons with 
experiments more extensively in Sect. VI. 

This paper is organized as follows: in Sect. II we illus- 
trate the method and the formalism; in Sect. Ill we pro- 
vide the normalization for the projected particle and hole 
excitations and in Sec. IV the probability for the particle 
and hole tunnelling; in Sect. V we check the accuracy of 
the GA using VMC computations. In Sect. VI we discuss 
the physical meaning of our findings and then we come 
to the Conclusions. In Appendix 1 we illustrate more 
specifically how to obtain the GA hopping factors intro- 
duced in Sect. II and present the full derivation of the 
analytical results of Sect. Ill and IV; in Appendix 2 we 
present some additional VMC results for projected BCS 
superconductors, in order to understand better some is- 
sues discussed in Ref. [l|. 



II. PROJECTED WAVE FUNCTIONS 

In this Section we will recall the crucial steps of the 
generalized scheme presented in Rcfs. [l|H and show how 
to adapt the formalism for the AFM wave functions. 

We consider an electronic system on a L-site lattice, 
where the position of each atom is given by Via — Ri+<5 Q : 
Rj belongs to the Bravais lattice and 5 a fixes the atomic 
site in the unit cell and then determines its sublattice. 
We will often use the shorthand notation Yi a = i + a and 
set the lattice parameter a=l. 

Considering states with a single particle (hole) excita- 
tion with respect to the ground state, it is convenient to 
work with the following partially projected wave func- 
tion, 
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Double occupancy is projected out for all the sites, but 
for the reservoir site l^f. In order to study AFM wave 
functions, the projector in Eq. ([2]) has the form 



whereas the parent fully projected state \*&) is 
|*)=P|*o), 

with the projector 
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where n; laa is the projected spin-cr electronic density for 
a site i on the sublattice a. The lattice is bipartite and 
divided into two sublattices, A and B: we fix as sublattice 
A the one with n a ^ > n Q j_. 

In order to work in the momentum space, we adopt the 
following convention to perform the Fourier Transform on 
a poliatomic unit cell: 



1 V p ikRi 



i(k+K)-S a t 



-k+K,cr' 



kK 



with k belonging to the Reduced Brillouin Zone and K 
to the reciprocal lattice. We can use a canonical trans- 
formation <j) to introduce the quasiparticle operators c*k 
from the original fermionic Ck operators: 



a lxa = X^ K (k<7)c k+K , 
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Then we can present Eq. ([6]) in an explicit form, using 
this Bogolubov transformation for the AFM: 



l kl<7 
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being J£af the AFM ordering wave vector. We con- 
sider the unprojected state |^o) which displays long- 
range AFM order, as found in previous studies: 21 ' 29 " — 



l*o> = IKjo) 
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|0) is the electron vacuum and 
1 , 
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where e k(T is the bare electronic energy spectrum and A 



(2) the staggering field. The condition ut 



1 comes 



from the normalization of the wave function ^o- 



Gutzwiller approximation 

The general form of the hopping renormalization fac- 
tors gtapa for a SDW (Spin Density Wave) is 



gtafia 
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where n° acr {n°^ a ) is the unprojected spin-cr density for a 
site i (j) on the sublattice a (/3) and n the charge density 
(the filling). Eq. (fTU| can be obtained from the general 
expression of the Gutzwiller factor g t , as we derive in 
Appendix 1A. If we consider the case of sites ia and j/3 
on different sublattices in an AFM, we obtain the factor 



1 



9tAB 
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2n^n Jn 



which gives the renormalization for the nearest-neighbour 
hopping (gtAB = 9tABa \ ^2 and are taken on the same 
sublattice). 

In order to calculate the excitation norms and the tun- 
nelling matrix elements, we work with expectation values 
such as 



(*; 7 i*; 7 > 
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that generalize the Gutzwiller renormalization scheme to 
partially projected wave functions. 

Thus we need the normalization (vf^ ) to evaluate 
the parameters g' 7 in Eq. To this aim, we introduce 

the quantity: 



X^ — 



Cd, T= o + C, 
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where Cd^=o is the number of states with empty or singly 
occupied reservoir site and Cd ly =i the ones with doubly 
occupied reservoir site. Then, X 7 is the ratio between 
the norm of the fully projected state and the norm of the 
state with one reservoir site. For the AFM the norms 
and (^ 7 |*[ 7 ) take the following form: 



(*i*>= E 

N at ,N al 

(K\%)= E 

N af ,N ai 
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JVof, N ail N h J \N t - N at , JVj. - N ah L/2-N + NJ a 
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JV«t, N ai , N ha J V^t " JV«r. N i ~ N <*1, L/2-N + N a 



L/2 



L/2 - 1 



N at , N ail N ha J \N t - N at -1,N±- N ai — 1, L/2 — N + N a + 1J a 



(13) 



where w a — e ^ m ( m - m °) ; being (3 a constant. w a is a sta- 
tistical weight associated to any electronic configuration: 
it is 1 for the homogeneous case and in first approxima- 
tion for the AFM in case of moderate magnetizations. 

The L-site lattice with N electrons (N a per spin a) is 
bipartite in two sublattices with L/2 sites. We indicate 
with N a the number of electrons for each sublattice (N aa 
per spin a) and with Nh a the corresponding number of 
holes, being Nf la = L/2 — ^}2 a N acr . For a given config- 
uration of the electrons on the two sublattices the total 
magnetization is M a = N a ^ — N a ^ — (L/2)m a with m a 
the correspondent magnetization density for that config- 
uration. Clearly, the possible values of M a are in general 
different from the thcrmodynamical total magnetization 
M = (L/2)m. In Eq. (fTB"]) the sums have the constraints 
M < N at <N t ;0< N ai <N i - M, with N a = N/2. 

For the homogeneus state (m=0), our Eq. (fT3"|) recov- 
ers Eq. (10) of Ref. H| in the thermodynamical limit. For 
the Neel state (m = n), all the spins f are on the same 
sublattice and all the spins J. on the other one: this im- 
plies that the double occupancy in the reservoir site is 
not possible and then = 1: in this case, the general- 
ized GA reduces to the standard GA and the auxiliary 
projector PL to the projector P. 



Eq. (|13p cannot be presented in a simple form, neither 
in the thermodynamical limit, and thus we check if a for- 
mula guessed after the homogeneous result can represent 
a reasonable approximation. The candidate form is 



(1 - n ht )(l - n hi ) 

which fulfils the known behaviour of A 7 both for m—0 
and m = n. Thus we investigate if Eq. (1141) approximates 
Eq. (|13j) adequately for intermediate m. 

At first, in Fig. [1] we perform an analysis of the finite- 
size effects affecting A 7 computed using Eq. (|13[) . for sev- 
eral fillings and for several lattices sizes (from L = 40 up 
to L ~ 10 3 ). Clearly the deviation of the obtained values 
with respect to the average ones reduces more and more 
upon increasing L and we conclude that X~ t computed 
on L ~ 10 3 site-lattices is representative of the result for 
L — > oo and will be compared with X® p evaluated using 
Eq. IH. 

Then, in Figs. [2] and [3J we compare systematically 
the behaviour of A 7 and X^ p . In Fig. [2] we consider the 
trends of X 1 as a function of m for several fillings and in 
Fig. [3] as a function of n for several magnetizations. The 
general outcome is that X® p is a very good approximation 
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FIG. 1: Factor Jf 7 as a function of 1/L, for the magnetization 
m — n/2 and several fillings n (L is the number of the lattice 
sites). 
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FIG. 2: Comparison of the factors X 1 (in symbols) and 
(in solid lines) as a function of the magnetization m, for sev- 
eral fillings n on a 1000-site lattice. 




FIG. 3: Comparison of the factors X 7 (in symbols) and X^ p 
(in solid lines) as a function of the filling n, for three magne- 
tizations m on a 4000-site lattice. 



of X 1 for small m and good for large m. In fact, the 
approximation is quite fair, because even for intermediate 
m (where the accord is less satisfactory), the quantitative 
discrepancy with the exact results is moderate (about up 
to 20%) and the qualitative trend is similar. 

We conclude that Eq. (fT4| is a controlled approxima- 
tion of Eq. (fT3|) for large L and therefore we decide to 
adopt Eq. (fT4")) for X 1 in the analytical derivations of the 
following Sections. 



Exact relations in the GA 

For the occupancy of the reservoir site Z7, the following 
exact expressions hold 

((1 -rcJrfX 1 - n hl))*'^ = X y (l-n h ) 
(n/ 7CT (l — n; 7) _ CT ))^j = X 7 ni 7cr 



= 1-X, 



(15) 



where (...)^ ee (^ 7 |...|^ 7 }/<^ 7 |^ 7 >^ 

In the following Sections we will use ci 7 as shorthand 
notation for (d; 7 )^j = ci; 7 . 



III. SINGLE-PARTICLE EXCITATIONS OF THE 
PROJECTED ANTIFERROMAGNET 



We consider the particle excitation 



\*t+K,<r) 

and the hole excitation 



^4+K,J*0> 



l*k+A>> = ^C k+K , CT |* ) 
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For the computations with 
mine the norms N^ +K a 



>, we need to deter- 
(*t + KJ*i +K J- We will 



k+K,a 



present here the results of our computations for the in- 
vestigated matrix elements where we used the GA; the 
full analytical derivation can be found in Appendix IB 
(particle excitations) and 1C (hole excitations). 



A. Particle excitation 

We intend to evaluate the norm of the projected par- 
ticle excitations: 



N, 



PPcUkJ^o)- (18) 



k+K,a 

After the due computations in the GA, Eq. (fT8"|) takes 
this form: 

AT+ -1 

= (1 - n) + -{gtA*n Aa + g t A-an A _ a ) 



N G 



9tAB(c k+K CT C k+K )* , 



(19) 
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FIG. 4: Ratio between the particle excitation norm JVjjJj. and 
the norm Ng of the correlated state as a function of the 
filling n, for several magnetizations m. 



where N G = = gtAAa [see Eq. JTUJ)] and n° Aa 

is the unprojected spin-er density for the A-sublattice. 
Thanks to Eq. |J7]), we can present Eq. (p^9|) in the follow- 
ing way: 



iV; 



1 



N G 



= (l-n) + -{gtAan Aa + g tA , 



r) 



N 



k+KAF,a 

N G 



9tAB(u k n klt7 + v k n k2a ) 
(!-«) + \{9tA- a n A _ c 

StAB(«k< lo . + u£< 2o .) 



9tAan Aa ) 

(20) 



kAcr 



kAo .a kAo .)ip , i.e. the unprojected mo- 
mentum distribution function associated to the magnetic 
quasiparticle state |^^ CT ) = o^a^IO)- 

For example, in Fig. |||we report the behaviour of the 
ratio between the particle excitation norm A+ and the 
norm N G of the correlated state as a function of the 
electronic filling n for some relative magnetizations m/n. 
This study helps us determine for which parameters the 
addition of one particle to the ground state is more or less 
favoured or even forbidden. Let us consider the case of 
an electron with k < kp for m = 0, being /cp the Fermi 
momentum (we will consider the same case for Figs. [5] 
and [5] as well). In such paramagnetic case, the particle 
excitation is not allowed for any fillings. Upon increas- 
ing the ratio m/n, this scenario is removed and particle 
excitations are possible for any finite doping. We can 
notice a typical trend associated to the increase of m/n: 
the probability is maximal for dilute systems and then 
decreases monotonically for larger fillings. The extreme 
situation occurs for the Neel state: the ratio N^/N G is 
constantly 1/2 until the doping closest to zero; for n = 1 
the probability for the particle excitation drops to zero 
and this is an isolated point. 

In general, the size of the magnetization affects a lot 
the entity of this excitation norm, but one feature persists 
for any value of m/n: it is forbbiden to have a projected 
particle-excited state at half-filling. 



J 



B. Hole excitation 



Here we determine the norm of the projected hole excitations: 



N, 



k+K,a 



N G 

After the due computations in the GA, we obtain 



(ttol4 + K, g -P^k+K,J*Q) 



(21) 
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N G 
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n m s,- a 1 — Xg 
~X~8 



(23) 



where we remark that I (to) belongs to the Bravais lattice and 7 (8) specifies the position in the unit cell and then 
the sublattice (see Sect. II). Finally, we find 
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(24) 
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FIG. 5: Ratio between the hole excitation norm and the norm Ng of the correlated state ^P, as a function of the filling n, 
for several magnetizations m. 



where the Z-factors are defined in Eq. (|2"51) . 



For example, in Fig. [5] we report the behaviour of the ra- 
tio between the hole excitation norm and the norm 
Nq of the correlated state as a function of the elec- 
tronic filling n for some relative magnetizations m/n. So 
we learn when projected hole excited states are possible 
and how the removal of one electron from the ground 
state is favoured or hindered by the change of the pos- 
sible parameters n and m. If compared to our previ- 
ous findings with N^/Nq, we can notice that the be- 
haviour of Nj^/No is only moderately affected by the 
choice of the parameters, except close to half-filling. In 
fact, upon increasing the magnetization, the probability 
for a hole excitation decreases moderately. At half-filling 
all the curves for different m/n diverge, as we could ex- 
pect thanks to the functional form of Eq. (PHI) and the 
knowledge that the A-factor is 1 at half-filling for any rel- 
ative magnetization. The trends revealed in Fig. 2] and 
[5] are opposite upon increasing m/n: the probability for 
projected particle-excited states is enhanced and the one 
for projected hole-excited states is reduced. 



Since 

N G N G 

= (*0|c k+ K^4 + K,J*0) jV*W 

N G N g ' 

we find that 

N + 

K-k* = ■ (26) 

Eq. (f2"6f implies that the particle tunnelling probability 
has the same properties described for the particle exci- 
tation in Sect. Ill A: it is enhanced for any fillings upon 
increasing the magnetization, but it is always forbidden 
at half-filling for any magnetizations. 



B. Hole tunnelling 



IV. TUNNELLING PROBABILITY 

A. Particle tunnelling 

Now we consider the tunnelling probability of a single 
particle into the projected state: 



Here we evaluate the tunnelling probability of a single 
hole into the projected state: 
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(* fc+if , CT |ck+K, ( x|*) 
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Nu + k^g 



(27) 
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After the analytical derivation reported in Appendix ID, 
(25) we find 
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^7<5<7 

For the AFM, Eq. (gHJ) becomes 
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where the Z-factors are defined in Eq. (f2"3")l and the Y-factors in Eq. 



(28) 



(29) 



(30) 



For example, in Fig.|6]we report the behaviour of the hole 
tunnelling probability Mj^. as a function of the electronic 
filling n for different relative magnetizations m/n. This 
figure illustrates the behaviour of the tunnelling probabil- 
ity of one electron extracted from the ground state, upon 
changing the filling n for some values m/n. In general we 
notice a decrease of the hole tunnelling probability, more 
pronounced for the Neel state. At half-filling both the 
hole tunnelling and the particle tunnelling are forbidden. 
This property of the half-filled AFM was expected for 




n 



FIG. 6: Hole tunnelling probability Af k ^ as a function of the 
filling n, for several magnetizations m. 



the introduction of a particle into the system and occurs 
also for his removal, displaying a particle-hole symmetry 
Upon increasing the magnetization, there is an enhance- 
ment of the tunnelling probabilities on the particle side 
and a reduction on the hole side with respect to the ho- 
mogeneous case presented in Refs. GUI. 

The particle and hole matrix elements presented in 
Eqs. (|26|) and (|30l) contribute to the conductance probed 
in the tunnelling experiments^ this illustrates the in- 
terest to work with the plane wave excitations |^ fe , K cr ). 
In fact, the excitations in Eqs. (fTS]) and (fT7| are respon- 
sible for the processes in tunnelling experiments, where 
the injected/extracted electrons are travelling waves. 

However, one can also be interested to perform a sim- 
ilar study, considering the projected Bogolubov excita- 
tions, i.e. the quasiparticle excitation 

\Hx*) = PaL,\* ) , (31) 
and the quasihole excitation 

l*fcA.) = P^xM- (32) 

The excitations in Eqs. (|3Tj) and (|32[) are associated to 
the electronic eigenstates and occur in photoemission ex- 
periments. 
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For the excitations [&^\ a ), one can evaluate the norms The latter quantities can be presented using the results 
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k\<T 



and the tunnelling probabilities 



formerly obtained for the excitations \^ 
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M k2a = 



*oKc] CCT Pe kCT P + <4 

C k+K AF ,cr 



(34) 



V. MONTE CARLO TEST 

In order to generate the possible electronic configura- 
tions for the AFM states, we proceed with the following 
scheme in our VMC computations. 

We use the uncorrelated state \^ ) obtained as the 
ground state of the mean-field Hamiltonian H AF : 



H AF 



E 

iajfjcr 



Uajf3c] a(T Cjp a . 



(35) 



where S? a is the z-component of the spin operator S on 
the lattice site ia. 

Then we consider the fully projected state \^f) and the 
partially projected state |^J 7 ) [see Eqs. from ((2]) to (0] 
on a two-dimensional square lattice, with the nearest- 
neighbour hopping ti a jfs = —t. In this case the bare 
dispersion is = — 2<(cosfc;i; + cosk y ) and the vector 

Kaf = (7T,7r). 

In our VMC computations, we evaluated the double oc- 
cupancy dry of the unprojected reservoir site for a compar- 
ison with the corresponding GA results: this will provide 
a test for the accuracy of the GA for projected AFM wave 
functions. In our Monte Carlo sampling of the electronic 
configurations, two kinds of updates are allowed to gen- 
erate a new configuration from the previous one: either 
one-electron hopping or spin-flip between two nearest- 
neighbour sites. The contribution of the spin-flip pro- 
cesses is more relevant close to half-filling and necessary 
at half-filling to update the configurations when double 
occupancies are not allowed, as in fully projected states 
'J. Both the types of update processes do not change the 
condition Nf = N± for the whole lattice but they allow to 
change the staggered magnetization of the microscopical 
configuration. Furthermore, since we need to evaluate an 
intensive local quantity (<i 7 ), we update the reservoir site 



oftener than the projected ones, in order to accumulate 
enough statistics in a shorter time. The standard update 
procedure uses the equiprobability to select two sites and 
then the configuration update interests the reservoir site 
on average once every L steps. If L « 10 3 the thermal- 
ization towards the equilibrium value will require a very 
long time for a local quantity. Instead, if we select the 
reservoir site 10-20 times oftener, this will shorten the 
time to achieve equilibrium value but nevertheless the 
statistics accumulated for the other sites is enough to 
have an undistorted description also of the surrounding 
environment. In this way we obtain an unbiased value for 
the double occupancy of the unprojected site, preserving 
the local detailed balance condition. 
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FIG. 7: Ratio between the magnetization m and the filling 
n as a function of the hole doping x for two values of the 
staggering field A on a 882-site lattice (VMC results). 
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Our VMC computations were performed typically 
on lattices ranging from 882 to 2500 sites: the size- 
dependent effects are quite negligible for our investiga- 
tions with AFM wave functions. In addition, we collected 
enough statistics in the number of electronic configura- 
tions generated for the sampling, so that the uncertain- 
ties of the data points here reported are within 1%, and 
then within the symbol size in the plots. 

In Fig. [7] we show the behaviour of the relative magne- 
tization m/n, upon increasing the hole doping x = 1 — n: 
this trend changes for different A, whose values are given 
in unit of the hopping parameter t. In our VMC ap- 
proach, we obtain a state with a finite magnetization 
upon applying a finite external field A. Thanks to this 
kind of study we find the whole correspondence between 
our control parameter A and the order parameter m of 
the AFM states. In fact, we know that for A = we 
recover the paramagnetic Fermi sea (m = 0), and for 
A = oo we obtain the Neel state with m/n — 1, while 
in the plot we present the response of the system for 
two intermediate values of A. For small A, m/n drops 
quite immediately upon increasing the doping and then 
is few percents of the half-filled value for intermediate 
and large dopings. This trend takes places because a 
weak field A is effective to induce a strong magnetiza- 
tion only if the electronic density is high, otherwise the 
system does not tend to place electrons with opposite 
spin on either sublattices. For large A, the decrease of 
the m/n upon increasing the doping is quasi- linear for a 
wide range of dopings: in a strong field it is convenient 
to place some electrons according to the AFM order also 
when only few electrons are on the lattice. A final inter- 
esting issue emerges in this investigation: for any nonzero 
A in the low-density limit the relative magnetization does 
not vanish and tends to a finite value. This implies that 
even diluted systems, and possibly the lattice gases, show 
quite sizable magnetization for finite not necessarily huge 
A: actually, the electronic system has a memory of the 
underlying bipartite lattice even if the vacancies exceed 
highly the occupied sites. 



Then we consider the double occupancy d 7 of the reser- 
voir site and we compare the results obtained in the GA, 
using Eqs. (0 and (fT5]). and the VMC results. For 
m = 0, we recover the results investigated in Ref. Q] for 
the projected Fermi sea, with an excellent agreement be- 
tween the GA and the VMC for any dopings. 

In Fig. [5] we plot the double occupancy e? 7 , as a func- 
tion of the doping, for two magnetizations: m = n (the 
Neel state) and m = n/2. For m = n, the probability to 
have a double occupancy is zero: in a very strong stag- 
gering field there is no chance to have two electrons on 
the same site, even on the reservoir one, and this is a 
behaviour shown both by uncorrelated magnetic systems 
and correlated ones. For intermediate staggering fields 
A and m/n, we find that the GA results are in an ex- 
cellent accord with the numerical ones for intermediate 



and large dopings. Instead, the agreement is less satisfac- 
tory for small dopings and particularly poor very close to 
half-filling. In general, this is quite reasonable due to the 
nature of the GA, whose performances are not very good 
for dense systems. In fact, in the GA all the configura- 
tions with the same number of f- and 4,-electrons give the 
same contribution to the Hilbert space, thus completely 
neglecting their spatial locations on the lattice. Then the 
GA does its best for homogeneous states and in general 
for not too dense systems: not so close to half-filling, the 
mutual interactions among the electrons on the lattice 
are quite sporadic and their specifical positions are not 
so crucial. 

The comparison of our results in GA with the corre- 
sponding VMC ones confirms the accuracy of the approx- 
imation. This is in general satisfactory for intermediate 
and large dopings and at least reliable for the qualitative 
trends for low dopings, where in any case the discrepancy 
is within 20% with respect to the numerical values. The 
agreement is worse for intermediate magnetizations. 

A similar study was performed in Ref. [l] for the pro- 
jected Fermi sea and for projected BCS-superconducting 
states. Since some open questions were left for the pro- 
jected superconducting states, we decided to perform the 
computations in Appendix 2, in order to clarify and pos- 
sibly fix these issues. 



VI. DISCUSSION 

The determination of the tunnelling matrix elements 
presented in Sect. IV is relevant to study the results 
of one-particle spectroscopies, as STM. 

These matrix elements are very difficult to be evalu- 
ated in an accurate realistic way, both analytically and 
numerically. The GA allows us to have insight into the 
structure of M k , and the accuracy of the approximation 
is tested in Sect. [Vj where we compare analytical results 




X 



FIG. 8: Double occupancy d 7 of the unprojected reservoir site 
as a function of the hole doping a; on a 882-site lattice: GA and 
VMC results are compared for intermediate magnetization m 
and for the Neel state (n is the filling). 



10 



in the GA and exact numerical VMC results (see for ex- 
ample Fig. [3]). The discrepancy between the approxi- 
mated and the exact results is very tiny for moderate 
and large dopings. The worst agreement is found close 
to half-filling for intermediate magnetizations, where the 
deviation is up to 20%. On the whole, the reliability of 
the approximation is satisfactory, and even very good for 
most investigated cases. 

Previous papers evaluated tunnelling probabilities for 
homogeneous states, mainly to study the spectral proper- 
ties of high-temperature superconductors^ - — Thus we 
decided to explore the tunnelling for the AFM states, not 
investigated yet. 

Tunnelling experiments allow for a complete study of 
one-particle processes: both the particle and the hole side 
of the spectrum can be obtained in the same experimen- 
tal setup, just inverting the potential difference between 
the metallic tip and the sample. Unfortunately, this kind 
of method can not be used for insulators because the po- 
tential barrier would be too large to have an acceptably 
high signal-to-noise ratio. The absence of STM spectra 
for insulators hinders us to have experiments for a direct 
comparison with our results. In principle, we need exper- 
imental results for correlated collinear AFM materials, 
not necessarily for insulators. But, in practice, all the 
strongly correlated materials which exhibit this specif- 
ical magnetic order are insulators, as the lowly doped 
cuprates i 38 ' 39 Other materials with moderate electronic 
correlations are AFM metals, as the iron-pnictides, but 
in this case the magnetic order is a stripe-like SDW. 40 
Our findings of Sect. Ill and IV could possibly interest 
other novel correlated materials. 

Single photoemission experiments can not describe 
both the hole and the particle side of the spectrum: in 
fact, hole excitations occur in direct photoemission, when 
the material is hit by photons and emits electrons; in- 
stead particle excitations take place in inverse photoe- 
mission, when the material is irradiated by an electronic 
beam and emits photons 4^ 

An important feature in tunnelling spectra for corre- 
lated materials is the particle-hole asymmetry42r— In 
uncorrelated materials the spectral weight for the exci- 
tations is symmetrical in its particle and hole side: this 
implies that the probability to add or remove one electron 
is the same. On the contrary, in correlated materials a 
particle-hole asymmetry is detected and the hole-side of 
the spectrum dominates. This is due to the strongly cor- 
related environment which penalises the introduction of 
an extra electron into the system but not its removal. Ac- 
tually, the tunnelling spectra show only a small asymme- 
try for low-energy excitation, being the most asymmet- 
ric features for high-energy excitations. Theoretically, 
the asymmetry was investigated so far especially for cor- 
related superconductors. The mostly accepted idea is 
that low-energy coherent excitations are responsible for a 
particle-hole symmetry in the spectrum also in correlated 
systems, and then the whole feature of the asymmetry 
would be due to the high-energy incoherent excitations. 



This scenario is illustrated by Randeria et a/4& using pro- 
jected BCS wave functions and a sum rule analysis for 
models with coulombic correlations; similar conclusions 
are also the outcome of other dedicated studiesi 44 i 47 ' 50 i 51 
The complementary idea is that a mean-field treatment 
of the electronic excitations can not capture any incoher- 
ent contribution to the total spectrum and then it would 
miss the asymmetry at all, as found in the explicit com- 
putations with the projected Gutzwiller excited states in 
homogeneous systems; 1 - A different position is expressed 
by Anderson and On g 48 i 49 and this is relevant for our 
present study. They deny that the Gutzwiller scheme for 
the projected quasiparticle excitations is intrinsically a 
mean-field approach limited to the coherent part of the 
spectrum and not able to find any asymmetry. Then 
in their approach they formulate projected excitations 
and obtain the particle-hole asymmetry in the tunnelling. 
This achievement is somehow in contradiction with other 
systematical investigations as in Ref. [l|. 

The previous studies are mainly motivated by the in- 
terest in the superconducting phase of materials like the 
cuprates. Instead, we decided to investigate the results 
of similar computations for projected AFM states. 

We point out that our framework is valid for low- 
energy particle (hole) excitations, the region where in 
the superconductors the tunnelling spectra exhibit only 
a small asymmetry. 43 

In Sect. Ill and IV, we determined the general prop- 
erties for the tunnelling elements, computed using the 
AFM projected wave functions defined by Eqs. (0]), (O, 
([5]) and © on a bipartite lattice, without further speci- 
fications of the model to be used. 

The main outcome of our work is that the increase of 
the magnetization has opposite effects on the behaviour 
of the particle and hole excitations norms (as shown in 
Sect. Ill) and also of particle and hole tunnelling matrix 
elements (as shown in Sect. IV). In fact, upon increas- 
ing the magnetization from the paramagnetic state to 
the Neel state, we found that the particle tunnelling is 
enhanced and, on the contrary, the hole tunnelling is re- 
duced though never suppressed. These changes in the 
tunneling are monotonous upon varying the magnetiza- 
tion. We emphasize that these findings are general for 
all the range of possible fillings. Therefore, this issue is 
expected to be robust, so to allow us for the following 
prediction: if one compares the tunnelling spectra of cor- 
related materials with collinear AFM order, one should 
find that the ratio of the particle-side vs hole-side sig- 
nal has to increase with the magnetization of the system. 
Actually, other effects contribute to give the conductance 
and they can possibly hide or alter this trend shown by 
the tunnelling matrix elements^; however, for materials 
with quite different magnetizations, the effect should be 
clearly observable in the low-energy sector of the spectra. 

Finally, the general form of our results in Eqs. (|20[) . 
(|2"rJl) . (1501 and ([Ml) is not explicitly particle-hole sym- 
metrical. However, we identified relevant particle-hole 
symmetrical cases, as the homogeneous limit for any fill- 
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ings or the half-filled case for any magnetizations. For 
practical purposes, correlated materials exhibit AFM or- 
der close to half-filling and then for this range of fillings 
our results can just allow for a tiny deviation from the 
particle-hole symmetry. 

VII. CONCLUSIONS 

Our investigation provides the tunnelling matrix el- 
ements for projected AFM wave functions in the GA. 
Such matrix elements contribute to the signal intensity 
in the one-electron spectroscopies for collinear AFM- 
ordered materials. To our knowledge, ours is the first 
theoretical work dedicated to study single-particle (hole) 
tunnelling for projected AFM states. 

We also tested the accuracy of our analytical results in 
GA through a comparison with VMC numerical results, 
findind a generally good agreement for any fillings and 
magnetizations . 

We explored the behaviour of the particle (hole) ex- 
citation norms and tunnelling probabilities for the fill- 
ings up to half-filling and for the whole range of possible 
magnetizations. The main outcome is that an increas- 
ing magnetization enhances the particle tunnelling and 
reduces the hole tunnelling. Then, we also identified a 
particle-hole symmetrical behaviour of the tunnelling for 
some relevant cases. 

For future investigations, we are developing a natu- 
ral extension of our present study and of the one in 
Ref. [l], evaluating matrix elements for projected two- 
particle (hole) excitations. This further analysis will be 
helpful to better understand the classes of experiments 
which probe a response function of the system, as for ex- 
ample charge and spin susceptibility in optical and neu- 
tron scattering, optical conductivity and resistivity. 
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Appendix 1 
A. Gutzwiller factors for the AFM 

The general expression for the hopping factor is 2 . 

( C Lt C J/3t> 

ffto/9t - —i — 

( c iaf C J/3t)0 

[n iat (l - n jP )n m {\ - n^)] 1 / 2 

For a SDW we know that 7ij Q — Tijp — n.l a — ixj^ — n and 
we can also relate the projected spin density rii aa and the 
unprojected one, n° ar7 : n iacr = Cn° lac7 {l - n° Q _ CT ), where 
C is a constant which we obtain, noting that 

and summing the following two equations: 

niaa — Cn° aa (l — n ia _ lT ) 
n>ia-<7 — Cn ia _ a {\ — n° iaa ) 



n 

Then we have this relation for a SDW: 

n iaa = -Z q o n iaa(^- ~ n iu-o)- (37) 

-y _ Ln iao n ia-tj 

n 

We can replace the latter relation for ni atT in Eq. ([36]) and 
obtain the factor g t for a generic SDW, given in Eq. (|10[) . 

To formulate the Gutzwiller factors for an AFM, we 
can adopt this kind of shorthand notation: riAa = 

TLB-cr = n<j. 



In an AFM also the magnetization m undergoes a renormalization with respect to the unprojected m . Using 
Eq. 437]), we find 



"it - n ii = 



(m ) 2 



2n 



If the sites ia and j/3 are on different sublattices, one obtains the factor gtAB = gtABo 

1 — n 1 — n 



gtAB = = 



1 - 2n°n0/n n 2 - (m ) 2 

2n~ 
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If they are on the same sublattice, one obtains the factors gtAAo = gtAa and gtBBo = gtBa- 

l-n 1 - n l _ l-n 2 - 7i + m° 

9tAt = 1 -2n^nl/nT^^ = n 2 - (m°) 2 2-n- m° ' 

2n 

_ l-n 1 - 71° _ l-n 2-n-m° 

9tBt = l-27i°7l°/nT^5 = n 2 _ ( m 0)2 2 _ „ + TO o • 

2n" 



Note that g tAa = g t B-<?- 



B. Particle excitation 



(*0|c k+K ,<r^4+K,«rl*0> 



L 
1 

I 



t , y,mS 



E<W 7 (1 - n l7 a)(l - n i7i _ CT )P/ 7 |* ) + E e lk - (i - m) e^ k+K H^^(* |P c ^ c t n5(T P|^ 

Z7 I'y^mS 



Ng 
L 



E 

/ 7 



(tf|(l-n/ 7 )|*) : y> eik .(( m ) i(k I. 



(38) 



where Afc = (^l^). For the diagonal contribution in the last step we have used the relation 

(*o|-P(l - nj 7 )P|*o) = (*o|-P/ 7 (l - - n /74 .)P/ 7 |*o)- Then we invoke the GA for the off-diagonal term in 

Eq. ([3HD: 



N k+K,a n s 1 



(*o| E -9*7^ e " "' *' 



ik-(i-m)„i(k+K)-( 7 -5)J 



mdrr l"yty 



E 5 *^ i5c 



a ik.(i-m)_»(k+K).(7-*)J 



C m5<T C i 7( rl*0) 



(i, 7 )=(m,<5) 



(l-n) 



(*0| - Eft77s- n !V 



; 7 



Z 7 ,m<5 



For the AFM the latter can be formulated in an explicit form thanks to Eq. ([7]): 



N, 



k+K,a 



(1 - n) + \ (gtAan° Aa + g t A-<jn\ _ a ) - #tAB(4 +K)ff c k+K> /*, 
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C. Hole excitation 



N, 



k+K,a 
No 



2 7 7 

i r i - x 

^ + 7 E 



l 



NqL 



E ik-(!-m) i(k+K)-(7-5)/ v r, ip/ t _ p' Itf, \ 



X 



N G L 



E e ik - ( '" m) e i(k+KH ''- 5) (^ m7 ,li CTCm5CT |^ m75 ) " ^ihM^h^lh-y) 

Z 7 ,m<5 Ij 



where P{ mi g = Yiia^i-y msO- ~ n ia,t n ia,i) ■ The last term in the above equation corresponds to a hopping process 
between two reservoir sites. The generalized GA assumes that the matrix elements are proportional to the square 
roots of the corresponding densities, as given in Eq. (fl~5|) . Invoking the GA, we obtain: 



N, 



k+K,a 



N G 



h 
1 



1 - X 



1 

y/n m 5<r(l - n m Sa) 

„0 



L [ 2^ x Xs e e (%. c m s*)o 

Z 7 ,m<5 

\j Xy(l — HJt)\/ X^ni^o- + XjUi^^a-^l — X^ 



\J Xj(l — n m g)\/ XsrimSa- + yj Xsn m s- a \Jl — Xj 



E 

; 7 



i 



'Z 7 (T 



Xy(l — n; 7 )-\/A 7 rii 7(J + -y/ X 7 ni 7 ^_ (T ^l — A 7 



fc+A'.cr 



= n a + 



I E [ V^l + I E [(e*-(-)e«(^K).h^) (c t ffCmJo _ z lSo 

£7 Irn-yS 



(39) 



with 



J -y8(T 



y(i - ^7) + 



' «; 7 ,-cr 1 — A 7 



^Z 7 o 



^m.5,-<T 1 — Xj 

X S 



D. Hole tunnelling 



Here we evaluate the tunnelling probability of a single hole into the projected state: 



M, 



(*£"+X,J<*+K,«r|*> 



(40) 
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We use Eq. (|T5j) and find: 



(*ol4 +K ^Pc k+Ki , 



N G 
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Appendix 2 

E. Monte Carlo results for projected 
superconductors 

As aforementioned in Sect. |Vl we performed compu- 
tations similar to the ones reported in Ref. 0], with the 
aim to understand better some topics discussed in that 
paper. For the superconductors, the lattice is not subdi- 
vided in two sublattices as for the case of the AFM, thus 
we drop the labels a and j3 to indicate the lattice sites in 
the formalism of this Appendix. 

We use the uncorrelated state |\&o) obtained as the 
ground state of the Hamiltonian H^ cs : 



HbCS — E tij c io c jt7 ~ M E C L 
ija ia 

+ E A ^ c ]i +h - c -) • 



(43) 



We consider projected BCS-superconducting states on 
a square lattice, with the nearest-neighbour hopping 



tij = —t. Then, we define our partially projected wave 
functions as 



I*!) 



with 



2(l + C/V(^) 2 + A^ 



(44) 



and with ££ 



e kcr H> 



M) being 



-2t(cos k x + cos k y 
H the chemical potential. In this Appendix the reservoir 
site is indicated with /. 

We check the accuracy of the GA for projected super- 
conductors for two different cases. Upon applying a finite 
A [see Eq. (|4"3"])]. we obtain BCS superconducting states 
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FIG. 9: Double occupancy di of the unprojected reservoir site 
as a function of the hole doping x for the projected s-wave 
superconductors on a 882-site lattice: GA and VMC results 
are compared for the field A = 2 in the case of onsite (sw) 
and nearest-neighbour BCS interaction (sw ext). 

and in particular we investigate what happens for the 
case of s-wave and d-wave symmetry of the gap. 

As we did in Sect. V for the projected AFM states, we 
consider the behaviour of the double occupancy di of the 
reservoir site in VMC in order to test the limits and the 
reliability of the GA. 



s-wave BCS superconductors. In this case the su- 
perconducting gap is isotropic. In Fig. [SJ we show the 
corresponding behaviour of the double occupancy of the 
reservoir site upon increasing the doping. In particular, 
we consider two different types of s- waves: the one given 
by onsite interaction A (recovering the results of Ref. [l]) 
and the extended s-wave given by nearest-neighbour BCS 
interaction. 

The trends of the double occupancy for the two types 
of s- waves are very different: for the extended s-wave the 
VMC results are very well approximated by the analyt- 
ical ones for all the dopings. On the contrary, for the 
local s-wave the double occupancy is highly enhanced 
with respect to the GA and for this specifical value of A 
the curve is quasilinear with the convexity tendentially 
downwards. 

These results imply that the GA fails completely as a 
good method to represent the behaviour of the projected 
onsite s-wave superconductor for any doping and this was 
quite unexpected, in that the GA should work well for 
the wave functions which do not exhibit spatial inhomo- 
geneities and should be extremely accurate in the diluted 
limit, while this is not the case. Nevertheless, we can ex- 
plain the physics responsible for this peculiar trend not 
captured within the GA. In fact, the onsite BCS interac- 
tion favours the onset of double occupancies with respect 
to the case of the intersite interaction which places the 
electrons preferentially on distinct sites. This scenario is 
clear for uncorrelated wave functions but actually it can 




x 



FIG. 10: Double occupancy di of the unprojected reservoir 
site as a function of the hole doping x for the projected d- 
wave superconductors on a 2500-site lattice: GA and VMC 
results are compared for the field A = 2. 



be exported also in presence of our Gutzwillcr projector; 
the different behaviour between the VMC trend for the 
pure s-wave and the trends for the extended one (in VMC 
and in GA) becomes robust upon increasing the size of 
A. In the limit of infinite applied field A, the double oc- 
cupancy of the reservoir site for the local s-wave will be 
maximal. In particular, for the uncorrelated local s-wave 
all the sites are unprojected and exhibit the double oc- 
cupancy d = n/2; in the strongly correlated case di = 1 
for A = oo. The GA is not sensitive to this limit for the 
projected local s-wave, as we can see in Fig. [S] the qual- 
ity of the approximation is highly unsatisfactory already 
for a finite value of A, quite far from the asymptotical 
regime of A = oo. 



d-wave BCS superconductors. The superconduct- 
ing gap shows this symmetry in the momentum space: 
Ak = A(cosk x — cosky). 

In Fig. I10[ we report the corresponding behaviour of 
the double occupancy di upon increasing the doping. For 
this case, the results presented in Ref. HI were not univo- 
cal, in the sense that the quantity ak=o — ltfe=oA>/!=o is 
indeterminate and thus the authors show the behaviour 
of the double occupancy for two possible values of ak=o'. 
and large, i.e. larger than all the other coefficients ak 
(see their Fig. 3). In general, they presented all their 
other results for ak=o large [as in their Fig. 2]. Actually, 
their choice was an ansatz and the quantitative discrep- 
ancy in the two different cases were quite relevant. In our 
computations, we obtain results compatible with theirs 
with larger ak=o and this supports their choice. 

Then one more open issue was left in the former VMC 
investigation by Fukushima et at: in fact they noticed 
that the trend of the double occupancy close to half- 
filling was extremely more sensitive to the size effects if 
compared to what occurs for lower fillings. This means 
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that we are to perform computations raising the num- 
ber of the sites L to several thousands in order to have 
a meaningful comparison with the GA, which is valid in 
the thermodynamical limit. According to the GA, di has 
to be 1 for n = 1: in Ref. [l| they found d 7 between 0.6 
and 0.8 for L between 1000 and 2000 sites. A further in- 
crease in the size is practically very difficult because nu- 
merical computations become long-time demanding, and 
a too limited statistics can cause not ergodic results. For 
L = 2500, di overcomes 0.8. According to the size-scaling 
trend, we are quite confident that this quantity is slow- 



ing tending to 1, but we can not claim this as a sure 
finding. This is the only case where we have a result so 
dependent in size also for very large lattices: typically 
the results for L > 500 are quantitatively very similar to 
the ones for other larger L. Clearly, one would be highly 
interested to fix this issue with VMC computations: a 
thermodynamical value of di ^ 1 for the projected BCS 
states in d-wave at half-filling would imply that the GA 
result is not realized in this case, while it is confirmed 
by the VMC computations for the other projected wave 
functions tested so far at half-filling. 
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